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Abstract: 

A two-dimensional lattice model for non-interacting fermions in a magnetic field with half a flux 
quantum per plaquette and iV levels per site is considered. This is a model which exhibits the 
Integer Quantum Hall Effect (IQHE) in the presence of disorder. It presents an alternative to 
the continuous picture for the IQHE with Landau levels. The large iV limit can be solved: two 
Hall transitions appear and there is an interpolating behavior between the two Hall plateaux. 
Although this approach to the IQHE is different from the traditional one with Landau levels 
because of different symmetries (continuous for Landau levels and discrete here), some charac- 
teristic features are reproduced. For instance, the slope of the Hall conductivity is infinite at 
the transition points and the electronic states are delocalized only at the transitions. 
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We consider in two dimensions non-interacting lattice fermions in a homogeneous magnetic 
field with half a flux quantum per plaquette. This problem was originally studied some time 
ago by Fisher and Fradkin [1]. It was shown by these authors that the large scale limit is 
equivalent to a two-dimensional Dirac theory. In a recent article Ludwig et al. [2] extended this 
model by introducing a staggered chemical potential. In the large scale limit this parameter 
appears as the mass M of the Dirac fermions. Therefore, the purpose of this parameter is to 
create a gap of the effective Dirac theory between the particle and the hole band. Thus two 
parameters control this system of non-interacting fermions: the staggered chemical potential 
and the energy of the particle E. In order to observe the integer quantum Hall effect (IQHE) 
one must introduce disorder. In the traditional approach to the IQHE, which is based on a 
continuum model with Landau levels [3,4], this is achieved by using random fluctuations of the 
chemical potential. However, for the model under consideration it was argued by Ludwig et 
al. [2] that randomness of the chemical potential alone cannot create a non-vanishing density 
of states (DOS) for M = E = 0. Their argument is based on a perturbation theory w.r.t. 
disorder. Consequently, disorder in M would not lead to a physical Hall transition because of 
the absence of states at the transition point M = E = 0. In order to get a non-zero DOS for the 
physical Hall system Ludwig et al. introduced another random quantity, an additional potential 

V which corresponds to energy fluctuation. The need of this potential is somehow surprising 
because a random chemical potential was sufficient to get the physical Hall transition in the 
famous approach to the IQHE by Pruisken et al. [3,4]. It will be argued subsequentially that, 
in contrast to the claim by Ludwig et al., a random mass is also sufficient for a non-zero DOS 
in the lattice model of Refs. [1,2]. 

First of all, we notice that the result of Ludwig et al. is in contradiction to an earlier work 
[5] where it was shown that the DOS is non-zero at M = E = 0. This result is based on a 
rigorous proof. Therefore, it would be interesting to understand why the perturbative approach 
of Ludwig et al. does not reproduce the correct behavior of the DOS. 

We will explain in this article that disorder in the Dirac mass M can lead to fluctuations in 

V (around V = 0). This is due to the fact that the same average Green's function can be created 
by different types of external fluctuations. Using iV-level fermions we obtain the fluctuations in 

V from a suitable representation and a non-zero DOS in the large N limit. 

To introduce the details of the model the pure system of non-interacting fermions on a two- 
dimensional lattice is considered along the lines of Refs. [1,2]. This describes lattice fermions 
with nearest neighbor hopping rate t = 1 and next nearest neighbor hopping rate t'/4 in a 
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staggered potential l) x+y . If we identify fermions with the four corners of the unit cell the 
related Hamiltonian reads in Fourier representation 

H(k) = 

I [i l + e~ ikx t(1 - e~ iky )(l - e~ ikx ) 1 + e~ iky \ 

l + e tkx —fi l + e~' lky —T(l — e~ lky ){l — e lkx ) 

—T(l — e lky )(l — e lkx ) 1 + e lky \x _i _ e *fc* 

\ 1 + e lky r(l — e' lky )(l — e~ lkx ) —\— e -^ _^ / 

(1) 

with t = it' /A. Expansion of k = (n,n) + ap for small p vectors leads to the large scale 
approximation which breaks up the Hamiltonian (1) into two independent Dirac Hamiltonians 
H± = a ■ p + o~3(m ± T") with Pauli matrices aj. The lattice constant a is scaled out with 
T" = t'/a and m = \ija. The two Dirac theories describe particles with different masses m±T ; , 
respectively. For large scale properties like the Hall transition it is sufficient to consider only 
the light particle with M = m — T' . The Hamiltonians obey individually a discrete symmetry 
H± — > —uzH±(jz for zero mass [2,5]. This is an important observation because the effective 
theory for the Landau level system obeys a continuous symmetry [3,4]. That means the two 
approaches to the IQHE are fundamentally different. Nevertheless, we will see that certain 
physical properties are similar in both cases. 

The dispersion of H- is E = ±^/M 2 +p 2 . This is a Dirac theory with a particle and a hole 
band which are separated by a gap with width 2M. We can vary the number of particles/holes 
by varying the width of the gap. Suppose, the energy is in the lower (hole) band. Varying 
the gap means we add/remove particles and holes to/from the edges of the gap. This has no 
effect on the current because those states are not occupied. On the contrary, we can have the 
energy in the upper (particle) band. Varying the gap means that we add or remove the same 
number of particles and holes. Consequently, there is again no net current change because 
both contributions cancel each other. The situation is different if the energy is inside the gap: 
only the lower (hole) band is completely filled whereas the upper (particle) band is empty. 
Varying the gap means adding/removing holes to/from the system if the energy passes a band 
edge. This implies an additional current. Finally, holes and particles can be exchanged by the 
transformation M — > —M. This gives a particle (hole) current for M > (M < 0), respectively. 
To make this more explicit the Hall conductivity a xy is calculated as the response to an external 
static field q y 



a 



xy 



j x /E y = 1- f Y^Tr[a x (H - iu + E)~^(H - iu + E + q y <J y )^ r ]^ (2) 

y r,r' 



where H is either H_ or H + and Y^' is the sum normalized with the number of lattice sites. 
Using the Green's function G(E — iui) = (H — iuj + E)~ x we obtain for q y ~ 



~iV '/ Tr[a x G{E - iu;r,r')G{E - iu;r' ,r")a y G{E - iu;r" ,r)] — . 



For infinite cut-off of the /c-integration (i.e. infinite bands) we find in units of e 2 



/n [2] 
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e(|M|-|£?|), 



(3) 
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where M can be the light or the heavy mass. This result reflects correctly the qualitative 
interpretation of the Hall conductivity. The Hall conductivity for the original lattice fermion 
problem is the sum of the Hall conductivities from the light and the heavy mass, such that the 
total a xy has a jump from to 1. 

Now we introduce disorder through a random Dirac mass M. Disorder creates a DOS inside 
the gap of the pure system. This means that the bands are broadened such that their inner tails 
can overlap. The gap is closed if the fluctuations SM are larger than M [5]. That means we 
have a compact region in M around M = for which the DOS p(E) is non-zero. We cannot 
distinguish particle and hole contributions to the current as in the pure Dirac theory because 
the particle hole symmetry is spontaneously broken in this case [5]: due to the fluctuations in 
M new complicated (mixed) states are created inside the gap of the pure system. They lead 
to a new Hall conductivity which interpolates between the two Hall plateaux. Once there is a 
gap again (if \M\ is larger than the maximal fluctuations |5M|) then the Hall plateaux appear 
(Fig.l). In the following we will generalize the model of Ref.[2] by introducing iV levels for the 
fermion Hamiltonian H±. Then it will be briefly discussed that the large iV Hall conductivity 
describes indeed such a behavior. Here we can use results obtained in previous studies of the 
Dirac Hamiltonian H- [8,9]. 

As a generalization of H± to iV levels of fermions we introduce H aa ' = Hq 01 ' — 5M^ a ' ' a 3 
(a, a' = 1, 2, N) with Hq 01 = H±5 aot and for the Hermitean random matrix SM we assume 



That means only random fluctuations couple the iV different Dirac systems. Physical quantities 
like the DOS per site or the conductivities are proportional to the number of levels N. It is 
natural to normalize these quantities by N. 




(4) 
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The purpose of the N —> oo limit is to obtain in a simple way a compact distribution to the DOS 
p(E, M) = (1/2) G^(E + i0 + ) from 5M which fills the gap. This is a consequence of the 
well-known fact that N — > oo gives a semicircular density p(M) = p(E = 0, M) [6]. However, for 
a non-compact distribution (e.g., Gaussian) and finite N there is no gap at all; the non-compact 
fluctuations 5M extend the DOS over the whole real axis. The trick is now to consider the 
large N limit and 1/N expansion (a finite order of these terms does not create particles outside 
the N — > oo DOS) as a technical device in order to maintain a density p(M) with compact 
support. There are Hall transitions at the edges of this 'band' which was created by the mass 
fluctuations. This result indicates that the Hall plateaux are only due to the opening of a gap in 
p(E, M). If we take N < oo there should be still a qualitative change in a xy due to a cross-over 
analogous to the transition with the compact DOS because the DOS in the gap is exponentially 
small for \M\ > M c . Therefore, o~ xy reaches almost the plateau value whereas we find a smooth 
change between the two plateau values by varying M between — oo and oo. Nevertheless, in a 
real system it is more natural to truncate the large fluctuations. 

The Green's function can be expressed formally as a functional integral for non-interacting 
fermions [7] 

G™;,(z;r,r') = [(H -5Ma 3 + zao)- 1 ]^.yy = -i [ *"' >M ,*" M exp(-Si) d$ r d® r (5) 

J r 

with the action (sum convention for a) 

S 1 = is z [-($, (H + za )$) + 6M? a '(*?' ■ a 3 $?)] (6) 

r 

with s z = sign(Im^) and the field = (*" M , Xr,^)- The first component is Grassmann and 
the second complex. The complex component is added to normalize the functional integral in 
(5). Averaging with Gaussian distributed fluctuations yields 

S 2 = -is z ($, (H + zao)$) + |p J^^r ■ o^rf- (7) 

r 

Thus we have derived an effective field theory for $ which serves as a generating functional for 
the average Green's function. It is important to notice that not only 5M creates the fermion- 
fermion interaction in (7) but also other types of disorder. For instance, the interaction can 
also be created by a term which couples to a matrix field (fx = 1, 4 includes the complex and 
Grassmann components): 

y r 
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This implies that the distribution 5M can be transformed into another distribution with a 
new 'random variable' Q (which does not have a probability measure but some generalized 
distribution including Grassmann variables). In other words, we can write 



([(H - 5Ma 3 + za )-T a )sM = ([(H - Q + zcx^T^q (9) 

The distribution which belongs to (...)q was investigated in detail in [8,9]. Here we present just 
the result for leading order in N: (...)q = / ... exp(— NS(Q, P)) Y[ r dP r dQ r with complex fields 
Q r , P r and 

S(Q, P) = - ^T[TY 2 (Q r 3 ) 2 + Tr 2 (P r a 3 ) 2 ] 
9 r 

+ log det(H - 2Q + za ) - log det(H + 2iP + za ) + 0(N~ 1 ). (10) 

The number of levels iV appears in front of the action. Thus the effect of disorder for N — > oo 
can be evaluated in saddle point approximation. The saddle point equation reads 

5 1 

— [-Tr(Q(7 3 ) 2 + log det(# - 2Q + za )] = 0. (11) 
oLJ g 

A second saddle point equation appears from the variation of P by replacing Q — > — iP. As 
an ansatz we take a uniform saddle point solution Q = —iPq = — (1/2) [27700 + M s a^\. Then 
(11) leads to the conditions 77 = (77 + 00 — iE)gI, M s = —Mgl/(1 + gl) with the integral 
I = J[(M + M s ) 2 + (77 + uj - iE) 2 + k 2 ]~ 1 d 2 k/2n 2 . This means disorder shifts the frequency 
uj — > uj + 77 and the Dirac mass M — > M' = M + M s , where 77(M, uj) and M 8 (M, u) are solutions 
of the saddle point eqn. The sign of 7/ is fixed by the condition that 77 must be analytic in uj. 
This leads to sign(7/) = sign(o;). Furthermore, p(M) is proportional to 77. The Hall conductivity 
per fermion level reads 



AP_ r 1 duj_ 

° xy ~ 2^ J (u + ri-iE- i\M'\)(u + i] - iE + i\M'\) 2^ ' ' " 1 



and with the approximation that M' and 77 do not depend on uj we get 

a xy w -\[arctan((M' + E)/rj) + arctan((M' - E)/r/)]. (13) 

The sum of a xy for both masses are plotted in Fig.l for E = 0. An analogous calculation gives 
for the longitudinal conductivity per level 

0xx ~ 4^2 K/2 - arctan((r/ 2 + M /2 - E 2 )/2r/|E|)]. (14) 
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The transition between the Hall plateaux does not occur at M = 0, as suggested in Ref. [2], 
but at M = ±M C where r](M) vanishes. The distance of these transitions is small for weak 
disorder and E = 0: ~ exp(— 2n/g). Thus it seems to be difficult to resolve the transitions in 
a real or numerical experiment. The slope of the Hall conductivity a xy (M) is infinite at the 
transition points in agreement with the Hall transitions for the Landau levels [15]. In the pure 
limit (rj — > 0) the conductivity a xx (M) is 1/47T in units of e 2 /h outside and zero inside the gap. 
This is consistent with the behavior of the Hall conductivity in the pure limit (eqn. (3)). In the 
disordered system a xx (M) vanishes also on the Hall plateaux where rj = provided E 2 < M' 2 . 
a xx (M) vanishes always for E = 0. This reflects localization of the states which are created by 
disorder at the center of the gap. 

Up to now only properties which are obtained in the iV — > oo limit were considered. Of 
course this is not sufficient to evaluate space dependent properties like the localization length. 
The latter requires the fluctuations around the saddle point. This can be studied on the basis of 
Ref. [9], where a divergent correlation length was found at the transition points M = ±M C . For 
M 7^ ±M C all length scales are finite which implies localization away from the transition points. 
A detailed investigation of the localized states must include the evaluation of the localization 
length exponent v. This exponent is known for the continuous system with Landau levels from 
semiclassical approximations [10] and numerical studies as « 7/3 [11-14]. However, other values 
were also found experimentally as well as theoretically [4,13,15,16]. 

In conclusion, a discussion of the IQHE for lattice fermions with half a flux quantum per 
plaquette and N levels per site was presented. In the large iV limit two Hall transition are found 
with an interpolating behavior between two Hall plateaux in contrast to the single transition of 
Refs. [2-4] . Although this system is not from the same universality class as the one described in 
Refs.[3,4] for the Landau level system, some properties are similar like the infinite slope of the 
Hall conductivity at the transitions. 
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Figure Caption 

Fig.l: Hall conductivity a xy in units of e 2 /h as a function of the staggered chemical potential 
M for disorder strength g = 2.2. 
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